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ABSTRACT 

We present a Chandra observation of the only low redshift, z = 0.299, galaxy cluster to contain 
a highly luminous radio-quiet quasar, H1821+643. By simulating the quasar PSF, we subtract 
the quasar contribution from the cluster core and determine the physical properties of the 
cluster gas down to 3" (15kpc) from the point source. The temperature of the cluster gas 
decreases from 9.0±0.5keV down to 1.3±0.2keV in the centre, with a short central radiative 
cooling time of 1.0 ±0.1 Gyr, typical of a strong cool-core cluster. The X-ray morphology in 
the central lOOkpc shows extended spurs of emission from the core, a small radio cavity 
and a weak shock or cold front forming a semi-circular edge at ~ 15" radius. The quasar 
bolometric luminosity was estimated to be ~ 2 x 10 47 ergs -1 , requiring a mass accretion rate 
of ~ 40 M yr~ 1 , which corresponds to half the Eddington accretion rate. We explore possible 
accretion mechanisms for this object and determine that Bondi accretion, when boosted by 
Compton cooling of the accretion material, could provide a significant source of the fuel for 
this outburst. We consider H1821+643 in the context of a unified AGN accretion model and, 
by comparing H1821+643 with a sample of galaxy clusters, we show that the quasar has not 
significantly affected the large-scale cluster gas properties. 

Key words: X-rays: galaxies: clusters — galaxies: quasars: individual: H1821+643 — inter- 
galactic medium — cooling flows 



1 INTRODUCTION 

The central radiative cooling time of the Intracluster Medium 
(ICM) in nearby galaxy clusters can drop below 1 Gyr and with- 
out a compensating source of heat the gas would rapidly cool down 
to low temperatures (for a review see Fabi an|1994) . However, the 
large reservoir of cool gas implied by this cooling flow model ex- 
ceeds the amount of molecular gas observed (eg. |Edge & Frayer| 
|2003[ >, and the inferred star formation rates by an order of mag- 
nitude (Johnston e et al.|1987||Hicks & Mushotzky|2005[|Rafferty[ 
|et al.|2 006; O'Dea et al.|2008[ >. In addition, recent high resolution 
X-ray spectroscopy has been unable to find the emission signatures 
of gas cooling below about 1 keV at the extreme rates predicted by 
the cooling flow model in the absence of heating ( Peterso n" et al.| 
[2003) [KaaTtra et al.||2004{ |Peterson & Fabian||2006^ . A heating 
mechanism is therefore required to stabilise the cooling in cluster 
cores (for a review see IMcNamara & Nulsen|2007) . 

High resolution images of nearby cluster cores taken by the 
Chandra X-ray Observatory have revealed complex structures, 
such as cavities, shock fronts and ripples, which have been pro- 
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duced by the central AGN (eg. Boehringer et al.|1 993 ; Fabian et al. 
[20001 120031 120061 IMcNamara et al.||2000 2001; Schindler et al. 



|2001||Forman et al.|20 05 ; McNa mara & Nu lsen 20071. X-ray cav- 
ities, or bubbles, expanded into the ICM by the AGN provide a 
direct and relatively reliable means of measuring the mechanical 
energy injected by the SMBH (eg. |Jones & De Young 2005}. In 
this way, AGN heating has been shown to be energetically capa- 



ble of balancing the cooling losses in cluster cores (Btrzan et al. 
2004l|Rafferty et al.|2006T,|Dunn & Fabian|2006| [McNamara et ~ 
2006). Repeated cycles of AGN mechanical heating could support 



the cluster gas in a relatively stable state. 

The galaxy cluster surrounding the quasar H1821+643, at 
redshift z = 0.299, was discovered optically by Schneider et al. 
{1992) and later detected in X-rays by Hall et al.H1997) and|Sax 



|ton et al.] ( |1997| l. The cluster is optically rich with an Abell rich- 
ness class > 2 ( [Lacy et aL]|1992) . The quasar's radio luminosity 
at 5 GHz lies in the classification for 'radio-quiet' quasars (at only 
10 23,9 WHz~' sr _1 ) and its 151 MHz radio lum inosity, which is 
10 25 3 WHz" 1 sr" 1 ^Blundell & Rawlings|200l} , is at the bound- 
ary observed to separate FR I and FR II structures. A deep VLA 
radio image of H1821+643 revealed that, although classified as a 
radio-quiet quasar, HI 82 1+643 hosts a giant 300 kpc FR I radio 
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source JBlundell & Raw lings 2001). HI 82 1+643 has previously 
been observed several times with the Chandra X-ray Observatory 
but in all cases with the gratings in place (HETG Fan g et al.|2002[ 
LETG Mathur et al. 2003 1 which have prevented a detailed analysis 
of the cluster emission. |Fang et al7| ( |2002) ) were able obtain a spec- 
trum of the cluster from the zeroth-order observation and found a 
temperature of ~ lOkeV and a metal abundance of 0.3Z;?). 

In this work, we present new Chandra observations of 
HI 82 1+643 without the gratings which have allowed us to analyse 
the underlying ICM and the interactions with the powerful central 
quasar. By comparing our results with a sample of nearby cool core 
clusters, we have explored the implications of two different modes 
of AGN feedback for HI 82 1+643 and investigated the impact this 
would have on the evolution of the ICM. 

We assume Hq = 70kms~ 1 Mpc~' , Cl m = 0.3 and £2 A = 0.7, 
translating to a scale of 4.4kpc per arcsec at the redshift z = 0.299 
of H1821+643. All errors are la unless otherwise noted. 
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2 DATA PREPARATION 

HI 82 1+643 was observed by Chandra for a total of 90 ks split into 
four separate observations taken over nine days with the ACIS-S 
instrument (Observation IDs: 9398 34ks, 9845 25 ks, 9846 18ks 
and 9848 1 1 ks). The data were analysed using CIAO version 4. 1 .2 
with CALDB version 4.1.2 provided by the Chandra X-ray Center 
(CXC). The level 1 event files were reprocessed to apply the latest 
gain and charge transfer inefficiency correction and then filtered 
for bad grades to remove cosmic ray events. For this analysis, the 
improved background screening provided by VFAINT mode was 
not applied. In the pileup affected region some of the real X-ray 
events would be marked as background events and excluded by the 
VFAINT filter leaving a hole in the quasar surface brightness. To 
filter the datasets for periods affected by flares we examined the 
ACIS-S 1 chip, which did not contain any extended cluster emis- 
sion. The background lightcurve for this chip was filtered using the 
LC_CLEAN scripQ provided by M. Markevitch and the ~ 4ks af- 
fected by flares were then excluded from the ACIS-S3 level 2 event 
files. The final cleaned exposure time was 85 ks. 

As the four separate observations were taken so closely to- 
gether, with effectively identical chip positions and roll angles, we 
were able to reproject them to a common position (Obs ID: 9398) 
and combine them. An exposure-corrected image combining the 
four final cleaned event files is shown in Fig. [I] The exposure map 
assumed a monoenergetic distribution of source photons of 1 .5 keV, 
which is approximately the peak energy of the source. Assuming 
that all X-ray counts within 3"radius of the point source were from 
the quasar and all counts external to that were from the cluster, we 
estimate there to be ~ 100,000 cluster counts and ~ 15,000 quasar 
counts, although the latter will have been significantly reduced by 
pileup. 

2.1 Pileup 

Pileup occurs whenever two or more photons, arriving in the same 
detector region and within a single ACIS frame integration time, 
are detected as a single event (see eg. Davis 2001 1. This effect hard- 
ens the source spectrum, because photon energies sum to create a 
detected event of higher energy, and causes grade migration. All 



Figure 1. Exposure-corrected image of the ACIS-S3 chip in the 0.5-7.()keV 
energy band smoothed with a 2D Gaussian o~ = 1.5". The logarithmic scale 
bar has units photons cm -2 s _ 1 pixel -1 . The straight line running across 
the entire image through the centre of the quasar is the readout streak, an 
instrumental artifact produced by the bright quasar. 

events detected by ACIS are graded based on the shape of their 
charge cloud distributions in a 3 x 3 pixel island. This grade is used 
to distinguish between a real photon or a background event, such as 
a cosmic ray. However, pileup modifies the distribution of charge 
and therefore alters the event grade. 

Following the CXC analysis of the Chandra PSrQ the ex- 
tent of pileup in the observation was estimated using the ratio of 
"good" grades (grades 0,2,3,4,6) to "bad" grades (grades 1,5,7) de- 
fined with observations from the ASCA satellite. The bad/good ra- 
tio is shown as a function of radius in Fig.|2]for two energy bands. 
The ratio was determined in circular annuli which scale logarithmi- 
cally in width to produce approximately the same number of source 
counts in each radial bin. 

Most of the emission from the quasar comes from the 0.5 — 
2.0keV energy band which, when piled up, was detected as higher 
energy photons above 2keV. The 2.0 — 5.0keV energy band was 
therefore used to estimate the radial extent of pileup for this ob- 
servation. Fig. [2] shows that pileup rapidly increases inside 3" and 
dominates at 1". The rise in the bad/good ratio beyond ~ 50" is 
caused by the increasing importance of the particle background. To 
minimise the effects of pileup, the region inside 3" was excluded 
from our analysis. This choice was also later confirmed to be appro- 
priate by repeating the MARX simulation of the quasar with pileup 
(Appendix A2). 



3 IMAGING ANALYSIS 

In Fig. [3] we show a merged exposure-corrected image of 
HI 82 1+643 centred on the bright, piled up quasar, with superim- 
posed 1.4 GHz radio contours from Blundell & Rawlings ( 2001| >. 
The extended cluster emission surrounding the quasar is elongated 
in a north-west to south-east direction. A surface brightness edge 
is also visible around the cluster core at a radius of approximately 
15" from the quasar. The edge, possibly indicative of a shock in the 



1 See http://cxc.harvard.edu/contrib/maxim/acisbg/ 



See http://cxc.harvard.edu/cal/Hrma/psf/ 
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Figure 2. Ratio of ASCA bad grades (gi57) to ASCA good grades (go2346) 
for H1821+643 in two energy bands. 

cluster gas, is most sharply defined around the north-west of the 
cluster core. 

The unsharp-masked image (Fig.[4]centre) highlights the edge 
in the surface brightness which can be seen to run round the cluster 
core from the north-west to south-east. This does not appear to be a 
symmetrical feature about the quasar; the radius of the edge varies 
from 13" to 17". Although the X-ray emission is extended to the 
north and south, there does not seem to be a close association with 
the 1.4GHz radio emission. The fractional difference image (Fig. [4] 
right) shows a bright extended spur of emission to the south-east of 
the AGN containing a surface brightness depression which could 
indicate a cavity. This X-ray depression is coincident with the in- 
ner part of the radio lobes but there is no close correlation with 
the outer region of the radio emission. Projection effects may have 
prevented a clean interpretation of the relation between X-ray and 
radio emission. 

There are potentially three additional cavities in the X-ray 
emission (north-west, north-east and south-west of the quasar, Fig. 
[4] right). However, these are not coincident with any 1.4 GHz ra- 
dio emission and it is not clear how they could have reached their 
current locations. 

To more quantitatively investigate the interaction between the 
quasar and the cluster, we extracted the properties of the cluster gas 
as close in to the quasar as possible. 



4 SPECTRAL ANALYSIS 

4.1 Separation of Quasar and Cluster 

Determining the properties of the ICM in the centre of the galaxy 
cluster requires a reliable separation of the quasar and cluster emis- 
sion. Although the cluster emission completely dominates over the 
quasar emission outside ~ 6" from the centre, we are particularly 
interested in probing the gas cooling times as close as possible 
to the quasar. Naively adding an absorbed power-law model, with 
freely fitting parameters, to account for this quasar emission gave a 
power-law normalization that was an order of magnitude below the 



expected value. The contribution of the quasar should be compara- 
ble to the cluster emission in the innermost radial bins suggesting 
that the majority of the quasar spectrum was incorrectly interpreted 
in the spectral fitting as hot cluster emission. We therefore gener- 
ated a simulation of the quasar PSF to explicitly account for the 
emission spilt-over from the quasar onto the cluster. 

The Chandra PSF has a complex structure that can be ap- 
proximately divided into two sections: a core produced by quasi- 
specular X-rays reflecting from the mirror surface and wings gen- 
erated by diffracted X-rays scattering from high frequency surface 
roughness. The diffractive mirror scattering also makes the PSF 
wings energy dependent: the high energy PSF is broader than at 
low energy. A detailed analysis of the Chandra PSF produced by 
the high resolution mirror assembly (HRMA) can be found on the 
CXC website0 

We used the Chandra ray-tracing program ChaRT |Carter| 
et al.|2003] > and the MARX softwar^] version 4.4, which projects 
the rays onto the detector, to generate simulated observations of 
the quasar. This allowed an analysis of the spatial and energy de- 
pendence of the PSF without the complication of cluster emission. 
ChaRT takes as inputs the position of the point source on the chip, 
exposure time and the quasar spectrum. However, because the core 
of the HI 82 1+643 observation is piled up, the quasar spectrum 
could not be determined simply from a small, quasar-dominated re- 
gion on the chip. Instead, the quasar spectrum was extracted from 
the readout streak. 

The analysis of the readout streak spectrum, a ChaRT sim- 
ulation of the quasar and a validation of this method for the test 
object 3C 273 are detailed in the appendix. In summary, although 
ChaRT does not provide an exact prediction of the observed PSF, 
the steep decline in the quasar contribution with radius compared 
to the cluster emission mitigated the effect of underestimating the 
PSF wings beyond 10". To test the accuracy of the ChaRT simula- 
tions, we repeated our method for a Chandra observation of 3C 273 
(obs ID. 4879), excluding the region containing the jet. The simu- 
lation accurately predicted the observed 3C 273 quasar spectra that 
were extracted in annuli centred on the nucleus. We therefore pro- 
ceeded with the ChaRT simulation of HI 82 1+643 and analysed the 
uncertainties introduced by this method in section |5T| 

4.2 Projected Radial Profiles 

Radial profiles of gas properties, such as temperature and density, 
were generated to analyse the gas properties in the core of the 
galaxy cluster. Spectra were extracted in a series of concentric an- 
nuli centred on the emission peak and excluding the readout streak 
and point sources (Fig. |3J. The radial bins were chosen to ensure 
a minimum of 3000 counts in each extracted spectrum with the 
minimum radius determined from the pileup analysis to be 3". This 
criterion ensured enough counts in each spectrum to provide a good 
spectral fit and constraints on the cluster properties. Point sources 
were identified using the CIAO algorithm WAVDETECT, visually 
confirmed and excluded from the analysis using elliptical apertures 
where the radii were conservatively set to five times the measured 
width of the PSF. All spectra were analysed in the energy range 
0.5 — 7keV and grouped with a minimum of 50 counts per spectral 
bin. The background was subtracted using a spectrum extracted in 
a sector at large radii, 170 — 210". Response and ancillary response 

3 See http://cxc.harvard.edu/cal/Hnna/psf/ 

4 See http://space.mit.edu/CXC/MARX/ 
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Figure 3. Left: Exposure-corrected image in the 0.5-7.0keV energy band smoothed with a 2D Gaussian a = 1, 5". The logarithmic scale bar has units 
photons cm -2 s~ ! pixel -1 . Right: with superimposed VLA 1.4GHz radio contours from Blundell & Rawlings (20oTJ. 
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Figure 4. Images of the central 0.8 x 0.8' with superimposed VLA 1.4GHz radio contours from Blundell & Rawlings 1 20011. Left: Exposure-corrected image 
in the 0.5-7.0keV energy band smoothed with a 2D Gaussian o~ = 1.5". Centre: Unsharp-masked image, created by subtracting images smoothed by 2D 
Gaussians with o~ = 1 and 5" and dividing by the sum of the two images (RMS noise ~ 0.1). The central 3" region containing the majority of the quasar 
emission was masked out before smoothing. Right: Fractional difference of each pixel from the average at that radius (RMS noise ~ 0.2). The resulting image 
has been smoothed with a 2D Gaussian o~ = 1.5". The SE X-ray depression is marked by a blue circle. 



files were generated for each cluster spectrum, weighted according 
to the number of counts between 0.5 and 7keV. 

These projected spectra were then fitted in XSPEC version 
12.5.0 (Arnaud 19961 with an absorbed power-law model, to ac- 
count for the quasar PSF, and an absorbed thermal plasma emis- 
sion model PHABS(POWERLAW) + PHABS(MEKAL) iBalu cinska^l 
Church & McCammon|[T992l |Mewe et al.][T985l [T9861 |Kaastra| 
T992l|Liedahl et al.|1995]>. Abundances were measured assuming 
the abundance ratios of Anders & Grevesse ( 1989 1. The parameters 
describing the quasar absorbed power-law model were fixed to the 



values given in Table [T] These values were determined by fitting a 
PHABS(POWERLAW) model to spectra extracted from the appropri- 
ate regions in the ChaRT simulation. The variation in the values of 
these parameters with radius is not caused by an intrinsic variation 
in the emission from the source. Instead, this variation is produced 
by the difference in the effective area of the detector for a small 
fraction of the quasar PSF compared to the cluster emission (see 
appendix). 

For the cluster model component, the temperature, abundance 
and normalization parameters were all left free. The absorbing col- 
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umn density was fixed to «h = 4 x 10 20 cm~ 2 determined from 
spectral fitting in the outer annuli. This value was consiste nt within 
error with the Galactic value «h = 3.44 x 10 20 cm -2 I Kalberla 



et al. 



19921. 



2005> . The redshift was fixed to z = 0.299 (Schneid er et al. 



4.3 Deprojected Radial Profiles 

A spectrum extracted from the centre of the cluster on the plane 
of the sky corresponds to a summed cross-section with a range of 
spectral components from the core to the cluster outskirts. To deter- 
mine the properties of the cluster core these projected contributions 
from the outer cluster layers can be subtracted off the inner spectra 
by making an assumption about the line of sight extent and depro- 
jecting the emission. 

In this work we used two deprojection routines: a straightfor- 
ward spectral deprojection (DSDEPROJ; |Sanders & Fabian||2007| 
|Russell et al.|2008] l assuming only spherical symmetry and an up- 
dated version of the 'X-ray surface brightness deprojection' code, 
which uses the additional assumptions of hydrostatic equilibrium 
and a suitable mass model, and was first described in Fabian et al. 
< fT98T) (see also |White et al.|1997[|ArTen et al.|2001||Schmidt et al 
[2001) . 

In more detail, DSDEPROJ starts with the background- 
subtracted spectrum extracted from the outermost annulus and as- 
sumes it was emitted from part of a spherical shell. This spectrum is 
scaled by the volume projected onto the next innermost shell (geo- 
metric factors from Kriss et al. 1983 1 and subtracted from the spec- 
trum of that annulus. In this way the deprojection routine moves 
inwards subtracting the contribution of projected spectra from each 
successive annulus to produce a set of deprojected spectra. The in- 
put projected spectra were analysed in the energy range 0.5 — 7keV 
and grouped with a minimum of 50 counts per spectral bin as be- 
fore. The radial bins were chosen to ensure a minimum of 3000 
counts in each deprojected spectrum and so are larger than the ra- 
dial bins for the projected spectra. The resulting deprojected spectra 
were also fitted with an absorbed power-law model to account for 
the quasar PSF and an absorbed thermal plasma emission model. 
The parameters describing the quasar model were fixed to the val- 
ues determined from the ChaRT simulation, where the ChaRT spec- 
tra were also deprojected (Table[TJl. The cluster spectral model pa- 
rameters were set as detailed in section |4~2] The spectra and best- 
fitting models for the inner deprojected spectra with significant con- 
tributions from the quasar emission are shown in Fig. [6] 

The X-ray surface brightness deprojection code uses the ob- 
served quasar-subtracted surface brightness profile, an NFW mass 
model (Na varro et al.| 1997) and assumes hydrostatic equilibrium 
to predict the temperature profile of the gas. The predicted temper- 
ature profile is then rebinned to match the DSDEPROJ result and the 
X 2 difference between the observed and predicted deprojected tem- 
perature profiles is calculated. The mass model parameters of con- 
centration and scale radius are iterated over using a Markov Chain 
Monte Carlo method to determine the best-fitting values which 
give the minimum % 2 '■ The MCMC routine was constructed with 
a Metropolis-Hastings algorithm comparing four chains of length 
190,000 with a burn-in length of 10,000. The metallicity was fixed 
to 0.4 Z r, for this analysis. 

This analysis assumes that the detected X-rays from this 
source were produced either by the quasar nucleus or by 
bremsstrahlung and line emission from the ICM. However, there 
could be a contribution from X-ray synchrotron or inverse- 
Compton emission produced by the radio jets and lobes. The 
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Figure 6. Spectrum and best-fitting model (solid line) for the inner DS- 
DEPROJ deprojected spectra requiring a quasar emission model: deproj fl , 
5 — 9" (upper), and deprojj, 9 — 15" (lower). The model is a combination 
of a power-law model (dashed line) to account for the quasar emission and 
thermal plasma model (dotted line) for the cluster emission. 



1.7GHz MERLIN observation of HI 82 1+643 showed that the ra- 
dio jets were sited within a radius of ~ 2", therefore an analysis of 
the X-ray emission at radii above 3" should not contain significant 
X-ray synchrotron radiation. The contribution of inverse-Compton 
emission, generated by relativistic plasma in the radio source lobes 
upscattering the jet-produced photons or the cosmic microwave 
background, was difficult to estimate. However, the spectra from 
3" were well-described by the two component spectral fits incor- 
porating the quasar and cluster emission and there is unlikely to be 
a significant extra contribution from inverse-Compton emission. 



5 RESULTS 

5.1 Cluster surface brightness profile 

The total surface brightness profile calculated from the exposure- 
corrected image for H1821+643 is shown in Fig. [7] The radial bins 
are l"wide in the cluster centre and this was increased at larger 
radii where the background becomes more important. The simu- 
lated quasar PSF was subtracted from the total surface brightness 
profile to produce the cluster surface brightness profile. The flatten- 
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Table 1. Table of best-fitting quasar parameters for each region. Column 1) Region name 2) Inner and outer radii of annulus (") 3) Galactic column density 
( cm -2 ) 4) Photon index 5) XSPEC power-law normalization ( photons keV~ 1 cm -2 s~ 1 at 1 keV) 6) X 2 I number of degrees of freedom 
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14/21 



deproj„ 5-9 0.06±0.01 1.57±0.05 3.8±0.2xl0~ 3 101/95 
deprojj, 9-15 0.006+g;g^ 6 L4llgjg 5.8+°j x 10~ 4 84/95 
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Figure 5. Exposure-corrected images in the 0.5-7.0keV energy band smoothed with a 2D Gaussian (J = 1.5". Left: The regions for projected spectra are 
shown by the magenta dashed lines (referred to as 1-19 from the centre outwards) and those for deprojected spectra are shown by the white solid lines (referred 
to as a-i from the centre outwards). The readout streak, shown by the green box, and point sources were excluded. Right: Zoom in of central region showing 
the annuli affected by quasar emission (only half of each annulus is shown for clarity). 



ing of the cluster profile from 15 — 20" corresponds to the surface 
brightness edge seen in Fig. [3] In addition, a second break in the 
cluster surface brightness profile can be seen at ~ 8". The cluster 
surface brightness profile decreases in the innermost radial bin sug- 
gesting that the quasar component could have been overestimated. 

Therefore, before proceeding with any analysis of these fea- 
tures, it was important to determine the uncertainty in the clus- 
ter surface brightness introduced by the quasar subtraction. This 
was achieved by comparing the observation and ChaRT simula- 
tion of 3C 273 as an estimate for the accuracy of the simulation 
of H1821+643. Fig. [7] shows a comparison of the ChaRT simula- 
tion of HI 82 1+643 with the observation and simulation of 3C 273. 
The surface brightness profiles for 3C 273 were scaled by the ra- 
tio of the fluxes of the two quasars, determined from the respec- 



tive readout streak spectra in the 0.5 — 7.0keV energy band. The 
ChaRT simulation was not expected to provide an exact prediction 
of the observed PSF, and indeed for 3C 273 the simulation under- 
predicts the observation, particularly in the PSF wings outside 10". 
The ChaRT model for the PSF wings is known to underestimate 
the wing^Jand therefore this underestimation is likely to also be a 
problem for HI 82 1+643. 

However in HI 82 1+643 the underprediction outside 10" is 
mitigated by the steep decline in the quasar contribution with ra- 
dius compared to the cluster emission. For the region inside 10" 
where the subtraction is important, the simulation of 3C273 un- 
derpredicts the observation by up to ~ 10%. This offset could be 



See http://cxc.harvard.edu/chart/caveats.html 
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Figure 7. Total surface brightness profile in the energy range 0.5 — 7.0keV 
for HI 82 1+643 (black circles), ChaRT simulation of the quasar surface 
brightness profile (blue crosses) and the cluster surface brightness profile 
produced by subtracting the quasar simulation from the total (magenta tri- 
angles). The surface brightness profile for 3C 273 (green squares) and corre- 
sponding simulation (white squares) are shown for comparison. The 3C 273 
profiles have been normalized by the ratio of the 0.5 — 7.0keV fluxes as 
determined from the readout streaks. The magenta dashed lines show the 
effect of increasing and decreasing the quasar contribution by 10%. 



due to the ChaRT model or alternatively caused by an underestima- 
tion in the modelling of the readout streak. Therefore, we estimated 
that the simulated quasar profile for H1821+643 could be offset by 
±10% inside 10" radius. This offset did not significantly affect the 
quasar spectral parameters determined from the simulated spectra 
extracted in large radial bins, which were used to determine the 
cluster properties close to the nucleus (Fig. 1 19|>. However, for the 
cluster surface brightness profile, the offset resulted in an uncer- 
tainty of 25% in the innermost cluster profile radial bin, dropping 
to 5% by the third bin and then an insignificant contribution above 
~ 6" (Fig. [7}. Therefore the decrease in cluster surface brightness 
observed in the central radial bin was not significant given the dom- 
inating contribution of the quasar. 

The cluster surface brightness profile was also deprojected and 
evaluated in sectors (Fig. [8} to study the core structures. The width 
of the radial bins was increased from those of projected radial pro- 
file to ensure enough deprojected counts in each bin and produce a 
reliable deprojection. The quasar contribution was first subtracted 
from the total surface brightness profile. Then the resulting cluster 
surface brightness profile was deprojected by assuming spherical 
symmetry to subtract the contribution from the outer annuli off of 
the inner annuli. The region from 3 — 5" was excluded from the de- 
projection because the shallow gradient of the core surface bright- 
ness profile produced large errors in the values, in addition to the 
uncertainty from the quasar subtraction. 

The SW sector (magenta) contains a section of the surface 
brightness edge at ~ 15" but no extended emission from the core. 
Fig.[4](right) suggests there could be a 'ghost' cavity from 5 — 12" 
in this sector. By analysing VLA 1.4 GHz and MERLIN 1.7 GHz 
observations, Blundell & Rawlings (2001) concluded that the jet 
axis in H1821+643 could be precessing. The 'tipped over' appear- 
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Figure 8. Exposure-corrected image in the 0.5— 7.0keV energy band 
smoothed with a 2D Gaussian a = 1.5". The labelled sectors correspond 
to the area used to produce the surface brightness profiles in Fig. [9] The 
two surface brightness edges are marked with circles at 8" (solid) and 15" 
(dashed). 

ance of the N plume of radio emission in Fig. [3] is suggestive of 
this. The X-ray depressions to the NE and SW identified in Fig. [4] 
could have been inflated by the jet at an earlier time, detached and 
risen buoyantly as the jet precessed. 

However the deprojected surface brightness profile does not 
appear to support this interpretation. Although there is less emis- 
sion from this sector, the surface brightness profile inside 14" is 
smooth without the clear decrement expected from a cavity. It is 
therefore likely that the fractional difference image (Fig. [4] right) 
highlights the excess emission in the extended arms compared to 
the remaining emission at the same radius. However, cluster sub- 
structure and deviations from spherical symmetry in the cluster 
core will cause problems for the deprojection routine, particularly 
at radii inside the surface brightness edge. We cannot rule out the 
possibility that these features are 'ghost' cavities (X-ray cavities 
without any detected coincident radio emission). 

The SE (green) and W (blue) sectors both contain extended 
arms of emission from the core and for the W sector this terminates 
particularly abruptly at the surface brightness edge at 15" (dashed 
line, Fig.[8]and Fig.|9j. For these two sectors there is also evidence 
for an inner surface brightness edge at 8" (solid line, Fig.[8]and Fig. 
[5). Behind this inner edge, there could be a radio cavity in the SE 
sector. However, although the surface brightness profile is consis- 
tent with a decrease at this radius, there are only two radial bins in 
this region and they contain additional uncertainty associated with 
the quasar subtraction. Therefore, it is not clear whether this is a 
radio cavity. 

5.2 Projected Spectral Analysis 

Fig.[l0]shows the projected radial temperature and metallicity pro- 
files for H1821+643 (also Table |2|. The temperature profile shows 
a cool-core with a broad range in temperature, almost a factor of ten 
from the cluster outskirts down to 1.3 ± 0.2keV in the centre. The 
temperature starts to decline steeply at ~ 18" and then breaks again 
at ~ 7" for a sharp drop into the cluster centre. The metallicity is 
approximately constant at ~ 0.3 Z Q at large radii but increases in 



© 0000 RAS, MNRAS 000, 000-000 



8 H.R. Russell et al. 



Radius (kpc) 
100 200 



10-' 



io- 7 




CQ 



10 



pq a 



io- 1 



10 J 



IO" 1 



ft c 
ft 



io- 1 



io- 1 




500 

i i i ( 

+ All 

4- SW sector 

-f- W sector 

-f- SE sector 



w 




10 



100 



Radius (arcsec) 



Figure 9. Top: quasar-subtracted surface brightness profiles in the energy 
range 0.5 — 7.0keV for H1821+643 extracted in sectors (Fig. [8j and the 
total surface brightness profile shown for comparison (black circles). Bot- 
tom: quasar-subtracted deprojected surface brightness profiles. The solid 
and dashed lines correspond to the circles in Fig. [8] 
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Figure 10. Projected radial temperature (upper) and metallicity (lower) pro- 
files. The quasar contribution has been subtracted. The superimposed lines 
show the effect of varying the normalization of the quasar component by 
±10% (dashed line) and ±20% (solid line). 



the cluster core to around ~ 0.5ZQ. The sudden drop in metallic- 
ity seen inside 6" is caused by an underestimate of the Fe K a line 
blend at ~ 5keV (Figure|6j. This is discussed in section[63] 

We examined the systematic error in the projected radial pro- 
files that could be introduced if the quasar normalization was varied 
by ±10% and ±20% and the power-law index varied by ±0.05. 
Fig. 



10 shows that within a radius of 5" a variation of ±20% in 



normalization of the quasar component would add an uncertainty 
of ±0.5 keV to the cluster temperature value and ±0.1 Z t :, to the 
metallicity. However, beyond 5" the contribution of the quasar PSF 
declines rapidly. Therefore this additional error is important only 
for the central two radial bins. Variation in the quasar power-law 
index did not produce significant additional error in the cluster pa- 
rameters. 

We applied spatially resolved spectroscopy techniques to pro- 
duce maps of the projected temperature and pressure in the cluster 
core (Fig. 111. The central ~4x4' was divided into bins using 



The central 

the Contour Binning algorithm ( Sanders 2006 1, which follows sur- 
face brightness variations. Regions with a signal-to-noise ratio of 
32 (~ 1000 counts) were chosen, with the restriction that the length 
of the bins was at most two times their width. Background and 
quasar spectra from the ChaRT simulation were subtracted from 
the observed dataset and appropriate responses and ancillary re- 
sponses were generated. Each spectrum was fitted with an absorbed 
MEKAL model with the absorption fixed to the Galactic value and 
the metallicity fixed to O.4Z . The spectra were grouped to con- 
tain a minimum of 20 counts per spectral channel and fitted in the 



energy range 0.5 — 7keV. The fitting procedure minimised the y}- 
statistic. The errors were approximately ~ 10% in temperature and 
~ 5% for the emission measure. 

Shown in Fig.[TT]are the emission measure per unit area, tem- 
perature and 'pressure' maps. The 'pressure' map was produced by 
multiplying the emission measure per unit area and the temperature 
maps. The emission measure map shows the strongly peaked core 
surface brightness, which flattens in the very centre and the over- 
all elongated morphology of the cluster in the NW to SE direction. 
The temperature in the cluster core drops down to 2keV and there 
is a region of cooler gas to the E forming part of a swirl, which 
could be similar to the spiral patterns seen in Abell 2029 ( Clarke 
|et al.|2004J , Ab ell 2204 (|Sanders et al.|2005J and in a temperature 
map of Perseus (Fabian et al. 2 006} . 

The temperature map also shows an increase to the NW from 
15 — 20" coincident with the strongest section of the surface bright- 
ness edge (Fig.|9](. This temperature increase is also found in a pro- 
jected temperature profile extracted in the NW sector (Fig.[8]l, how- 
ever this increase is not statistically significant enough to draw any 
further conclusions on the nature of the surface brightness edge. In 
addition, the pressure map does not show a particularly sharp drop 
from 15 — 20" which would indicate a shock. 

It is not obvious from any of the maps shown in Fig. [TT] in 
which direction the main radio axis is oriented. Apart from a small 
low pressure region to the SE, which could relate to the small ra- 
dio bubble, the pressure map appears symmetric about the central 
source with no indication that it contains a large FR I source. 
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Figure 11. Left: emission measure per unit area map (units cm~ 5 arcsec~ 2 ). The emission measure is the XSPEC normalization of the MEKAL spectrum 
K = Elf (4 x \0 lA 7lD 2 A (1 + z) 2 ), where EI is the emission integral El = J n e nudV. Centre: temperature map (keV) with arrows marking the possible swirl of 
cool gas. Right: 'pressure' map (keVcm~ 5 arcsec~ 2 ) produced by multiplying the emission measure and temperature maps. The maps cover the same spatial 
scale and the quasar contribution has been subtracted. The readout streak and central 3" radius region affected by pileup have been masked out. 



Table 2. Table of best-fitting cluster parameters for each projected region. 
Column 1) Region name 2) Inner and outer radii of annulus (") 3) Temper- 
ature (keV) 4) Abundance (Zq) 5) % 2 I number of degrees of freedom 



Region 


Radius 


Temp. 


Abund. 


X 2 /doi 


proji 


3.4-4.2 


1.3 ±0.2 


I0+ 011 
u - lu -0.06 


15/24 


proj 2 


4.2^1.9 


2 0+ 03 


2+ 02 
""'-0.1 


26/25 


pr°j 3 


4.9-5.7 


9 7+0.4 
z - '-0.3 


06+° 16 
u - uo -0.06 


26/25 


proj 4 


5.7-6.4 


3.0±0.3 


5+ 03 


26/27 


pr°j5 


6.4-7.1 


3.5±0.3 




29/26 


prqje 


7.1-7.9 


o 4 +o.3 
J -^-0.2 


0.8i«;| 


27/27 


proj 7 


7.9-9.3 


3.9±0.2 


0.5±0.1 


38/51 


prqjs 


9.3-10.8 


4 Q+0.4 
+ y -0.3 


0.3±0.1 


57/50 


proj 9 


10.8-12.3 


c 7+O.6 
•'•'-0.5 


0.3±0.2 


22/41 


prqjio 


12.3-13.8 


« 7+0.6 
3 -'-0.5 


0.8±0.2 


34/36 


prqjll 


13.8-15.7 


7 9+0.8 
'• z -0.7 


0.3±0.1 


44/44 


proj,2 


15.7-17.7 


8±1 


06+° 16 
u - uo -0.06 


42/40 


projo 


17.7-19.7 


7 4+0.8 
'•^-0.7 


0.6±0.2 


38/39 


proj,4 


19.7-23.6 


o 9+0.7 
°- 2 -0.6 


0.4±0.1 


90/75 


projis 


23.6-29.5 


9.0±0.6 


0.24 ±0.09 


119/110 


projis 


29.5-39.4 


8.3±0.4 


0.30 ±0.07 


166/150 


projn 


39.4-59.0 


8.6±0.4 


0.22 ±0.06 


206/183 


projis 


59.0-98.4 


q 7+0.6 
y -'-0.5 


0.28 ±0.08 


236/214 


projig 


98.4-147.6 


12+3 


1+ 03 
u - 1 -0.1 


195/209 



5.3 Deprojected Radial Profiles 

Fig. [12] shows the deprojected radial profiles for both the sur- 
face brightness and spectral deprojection methods applied to 
H1821+643. The results for the spectral deprojection are also listed 
in Table [3] The surface brightness deprojection produced a mass 
profile with a concentration c = 1.3^22 and r^joo = 2.5^Q :)Mpc. 
The spatial region from 3 — 5" was excluded from the deprojection 
analysis because the shallow gradient of the core surface bright- 
ness profile and the core cluster substructure produced large errors 



Table 3. Table of best-fitting cluster parameters for each deprojected region. 
Column 1) Region name 2) Inner and outer radii of annulus (") 3) Tempera- 
ture (keV) 4) Electron density (cm~ 3 ) 5) Abundance (Zq) 6) x 2 I number 
of degrees of freedom 



Region Radius 



Temp. 



Abund. 



X 2 ldof 



(xl0~ 3 ) 



deproj fl 


4.9-8.9 


2.4±0.2 


46 ±2 


4 +0 - 2 


72/58 


deproj 6 


8.9-14.8 


4 3+ - 4 


26 +0 - 6 


0.2±0.1 


71/58 


deproj c 


14.8-24.6 


7 7 + 1.0 
'•'-0.8 


11.7±0.3 


0.5 ±0.2 


64/58 


deproj rf 


24.6-34.4 


7 7 + 1.0 

'•'-0.8 


8.8±0.2 


0.5 ±0.2 


66/58 


deproj f 


34.4-51.7 


7 6 +0 ' 8 


5.41 ±0.09 


0.4 ±0.2 


69/58 


deproj/ 


51.7-64.0 


10+ 3 


3.8±0.1 


2 +0 - 3 


51/58 


deproj s 


64.0-88.6 


Si s + '-O 


2.50±0.05 


0.3 ±0.2 


46/58 


deproj;, 


88.6-157.4 


12+ 2 


1A,S -0.03 


1 +0 ' 3 


45/58 



in the parameters, in addition to the uncertainty from the quasar 
subtraction (Fig. [7}. 

The deprojected temperature profile shows a steep drop in- 
side the central 100 kpc, from 10^ keV at ~ 250kpc down to 
2.5^Q2keV in the cluster centre. The surface brightness deprojec- 
tion result is similar but shows an increase in temperature within 
the central 40 kpc. This appears to correlate with a decrease in the 
electron density, calculated from the emission measure of the de- 
projected spectrum, however the errors are consistent with a flat 
profile. Although there is a break in the density profile at ~ 15" 
corresponding to the surface brightness edge, there does not appear 
to be a sharp discontinuity to indicate a shock. A jump in either 
temperature or density cannot be discerned for the DSDEPROJ re- 
sult because large radial bins were required to ensure enough counts 
for a good deprojection. The assumption of hydrostatic equilibrium 
in the surface brightness deprojection would also have smeared 
out any discontinuity associated with a shock. Finally, the surface 
brightness edge is not radially symmetric about the quasar. There- 
fore the circular annuli centred on the quasar that were used for 
the deprojection will have inevitably smeared the spectrum from 
this region over neighbouring annuli. Deprojecting the cluster in 
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the NW sector (Fig. [8j where the edge is strongest and narrow- 
est would be more likely to produce a conclusion about its nature, 
however there are not enough X-ray counts to deproject in sectors. 
It was therefore not possible to determine unequivocally the nature 
of the surface brightness edge with the available data. 

The metallicity profile produced by spectral deprojection is 
consistent with a constant value of ~ 0.4ZQ in the core, decreasing 
to 0.2 — 0.3 Z at large radii. The assumption of a constant abun- 
dance of 0.4Zp ; for the surface brightness deprojection was there- 
fore deemed to be suitable. 

Both deprojections were repeated with a ±10% variation in 
the subtracted quasar contribution. This resulted in an uncertainty 
of ±0. 1 keV for the innermost annuli in the surface brightness de- 
projection, smaller than the statistical error of ±0.3keV. The vari- 
ation produced a ±0.2 keV uncertainty for the innermost bin of the 
DSDEPROJ result but was insignificant at larger radii. The increase 
in the central temperature shown in the surface brightness depro- 
jection is therefore an inconsistency between the two deprojection 
methods. A finer spatial binning over the inner radii was tried for 
the spectral deprojection in order to confirm this feature. However, 
the flat surface brightness profile in the core resulted in a larger 
fraction of projected emission in each annulus. Therefore, increas- 
ing the number of annuli to be deprojected produced much larger 
uncertainties in the central shell. The increase in core temperature 
is not seen in the projected temperature profile (Fig, 1 10} or the tem- 
perature map (Fig.[TTJ. Whilst the lack of confirmation in projected 
profiles or the alternative, low resolution deprojection method does 
not preclude an increase in core temperature, the feature could in- 
stead be an artifact of the deprojection method. 

The additional assumptions of hydrostatic equilbrium and an 
NFW mass model made by the surface brightness deprojection 
method might not be appropriate for the cluster core. For a sam- 
ple of fourteen galaxy clusters observed with Chandra, |Voigt &| 
Fabian (2006) found that the mass distribution for the majority of 
the objects was consistent with the NFW model, however four ob- 
jects in the sample exhibited a flatter core. Forcing an NFW mass 
model in the cluster core when the gas mass was in fact much lower 
could, via hydrostatic equilibrium, artificially boost the calculated 
ICM temperature. Structural features in the cluster core, such as 
a shock or dense blob of gas, are also expected to cause spurious 
deprojection results. 

By combining the deprojected temperature T and electron 
density n e , we calculated derived properties such as the electron 
pressure P e = kgn e T, radiative gas cooling time 
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where A(T) is the cooling function, and entropy S = kBTn e 
We assumed the temperature and density were independent of each 
other in order to calculate errors on the derived quantities. The ma- 
jority of the error was in the temperature values so this was a rea- 
sonable assumption. 

Fig.[T3]shows the electron pressure increases steadily towards 
the cluster core. However, beyond the surface brightness edge the 
DSDEPROJ result flattens possibly to a constant pressure within 
50kpc. The surface brightness deprojection result assumes hydro- 
static equilibrium and an NFW mass model, which constrain the 
pressure profile. Therefore this conflicting result may indicate that 
the radio source provides a significant source of non-thermal pres- 
sure in the cluster core. The value of the pressure for the innermost 
radial bin in the spectral deprojection result should be considered 
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Figure 12. Quasar-subtracted deprojected temperature (top), electron den- 
sity (centre) and metallicity (bottom) profiles for the spectral deprojection 
method (red diamonds) and the surface brightness deprojection method 
(black circles). The radius extends in to 3"to enable a comparison with Fig. 
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a lower limit on the total pressure. It is therefore difficult to deter- 
mine whether there is a jump in pressure across the surface bright- 
ness edge at 13 — 17" indicating a shock. If instead the pressure 
smoothly changes across this region the feature could be a cold 
front (see Markevitch & Vikhlinin 2007 and references therein). 

The gas cooling time falls steadily down to l.Ojyj Gyr at 
a radius of 7"(30kpc) in the core of the cluster (Fig. 14 1. With 
such a short central gas cooling time, we would expect that gas 
might be cooling down below X-ray temperatures. The cooling ra- 
dius was defined to be the radius at which the cooling time fell to 
7.7Gyr, the time since z= 1. The cooling radius of H1821+643 is at 
20"(90kpc). The X-ray mass deposition rate inside the cooling ra- 
dius was calculated by using the XSPEC MKCFLOW model, which 
models the emission from gas cooling between two temperatures 
where the normalization is the mass of cooling gas. 

The DSDEPROJ deprojection process was therefore repeated 
using a larger innermost annulus 5 — 20" (22 — 90 kpc). How- 
ever, the majority of the mass deposition in clusters occurs in the 
cluster centre, where the radiative gas cooling times are short- 
est. The crucial inner 20 kpc was necessarily excluded from our 
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Figure 13. Quasar-subtracted deprojected electron pressure profile for the 
spectral deprojection method (red diamonds) and the surface brightness de- 
projection method (black circles). The shaded region denotes the radial ex- 
tent of the surface brightness edge from 13—17". 



analysis, therefore to get an upper limit on the mass deposi- 
tion rate, we assumed a single-phase cooling flow and fitted 
a PHABS(POW)+PHABS(MEKAL+MKCFLOW) model to the cool- 
ing region spectrum in XSPEC. The lower temperature of the 
MKCFLOW component was therefore fixed to 1 keV and the higher 
temperature was tied to the temperature parameter in the MEKAL 
component. The abundance was fixed to O.4Z . The quasar spec- 
tral properties were calculated as before by fitting to the ChaRT 
simulation of the cluster spectrum from the cooling region. The 
spectral fit resulted in an upper limit on the mass deposition rate 
of 300 ± 100M©yr within the cooling radius. Direct detections 
of the amount of gas cooling out of the X-ray requires deep high 
resolution spectroscopy (eg. Pete rson et al.|2003) . 

The heating rate refers to the spatial distribution in heat re- 
quired to balance the cooling of the cluster gas and was calcu- 
lated by dividing the total X-ray luminosity in each shell by the 
shell width. Fig. [14] shows that the rate at which energy must be 
deposited in shells of equal width is approximately constant at 
10 43 ergs -1 kpc -1 within lOOkpc, decreases in the radial bin as- 
sociated with the surface brightness edge, and then drops off be- 
yond 200 kpc. The total required heating rate to balance the cooling 
losses within the cooling radius is therefore ~ 10 45 ergs -1 . 



5.4 Potential cavity heating 

Assuming the X-ray depression at a distance of 24 kpc from the 
nucleus is a radio cavity (Fig. [4}, we follow the method of |Dunn &| 
Fabian (2004) (see also Birzan et al. 2004) to calculate the cavity 
heating power. The energy input of the cavity is the sum of the 
cavity's internal energy and the work done, 

1 y 

E = -PV + PdV m —t-rPV, (2) 

y— 1 y- 1 

where P is the thermal pressure of the surrounding ICM, V is the 
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Figure 14. Quasar-subtracted deprojected ICM cooling time (top), entropy 
(centre) and heating rate required per kpc width shell (bottom) profiles for 
the spectral deprojection method (red diamonds) and the surface brightness 
deprojection method (black circles). The radius extends in to 3"to enable a 
comparison with Fig. |10| 



volume of the cavity and y is the mean adiabatic index of the fluid 
in the cavity. The relativistic case of y = 4/3 was assumed. 

We used two different methods to place upper and lower lim- 
its on the potential cavity heating power. For a lower limit, we as- 
sumed there is no additional non-thermal pressure in the cluster 
core and used the DSDEPROJ weak pressure gradient shown in Fig. 
[T3| For an upper limit, the gas pressure profile from 40— 150kpc 
was extrapolated inwards using a powerlaw to produce an estimate 
of the total central pressure, similar to the surface brightness de- 
projection result. The instantaneous mechanical power for the cav- 
ity is /^ u bble — 4PV / // a g e , where ? age is the bubble age. The cavity 
was assumed to be spherical with a radius of ~ 9 kpc, estimated 
from the image using a circular aperture, and located at a distance 
of ~ 24 kpc from the cluster centre. Two estimates of the cavity age 
are commonly used: the time required for the cavity to rise to its 
current location at the speed of sound, f soun d, and the time required 
for the cavity to rise buoyantly at its terminal velocity, ?buoy The 
sound speed timescale is more applicable to young radio cavities 
recently inflated by the AGN. Its value was 0.03 Gyr, which gives 
a cavity heating power of 1 — 2 x 10 44 ergs -1 . 
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Although there could be extra AGN heating power produced 
by weak shocks and sound waves (Fa bian et al.|20 06 1, cavity heat- 
ing in HI 82 1+643 is not currently supplying enough heat to replace 
the majority of the ~ 10 45 ergs _1 radiated by the cluster gas. The 
possible 'ghost' cavities could potentially produce an additional 
power of 9 x 10 44 erg s~ 1 , calculated using the buoyancy timescale. 
However, it is not clear from the deprojected surface brightness pro- 
files (Fig. [9} that these features are cavities. 



6 DISCUSSION 

6.1 Evidence for quasar and cluster interactions 

The morphology in the central lOOkpc of HI 82 1+643 suggests a 
complex interaction between the quasar and the surrounding clus- 
ter gas. There are extended arms of emission and an inner radio 
cavity with a bright rim in the core and a weak shock or cold front 
surrounding most of the cluster core at ~ 15". However, the mech- 
anisms by which the quasar output is coupled to the ICM are far 
from clear. 

The extended arms of X-ray emission to the N and SE appear 
to be correlated with the main axis of the extended FR I radio emis- 
sion. These X-ray features may have been pushed or dragged out 
from the cluster core by the expanding FR I structure. Although the 
inner radio cavity is too close to the quasar to be analysed effec- 
tively, this indicates that there could be some mechanical feedback 
in the cluster core. The ends of the FR I structure appear to termi- 
nate at the surface brightness edge seen in the X-ray so this feature 
could be a weak shock generated by the radio source expanding 
into surrounding cluster gas. However, Blundell & Rawlings ( 2001 1 
found that the outermost ends of the extended radio emission do not 
appear to end in shocks. They inferred that the extended emission 
has not grown any faster than the ambient sound speed of the ICM. 
In addition, the extended X-ray structure to the W is not associated 
with any radio emission and the surface brightness edge appears 
strongest in this radial direction (Fig. [9j. 

Based on optical imaging, Hutchings & Neff ( 1991 ) suggested 
that HI 82 1+643 is in the late stages of a mild tidal event or merger, 
although there is no sign of a disturbing object. |Fried| {1998} 
found that the line widths of the extended emission-line gas in 
H1821+643 produced velocity dispersions in agreement with this 
interaction scenario. However, velocity dispersions of a few hun- 
dred kms~' are also consistent with the kinematics of emission- 
line nebulae in cool-core clusters ( Hatch et al.| |2~007l. The X-ray 
features could be interpreted as resulting from a merger interac- 
tion which has disturbed the gas in the cluster core. |Blundell"&| 
|Rawlings| ( [2001[ > estimated that if the radio plumes have emerged 
no faster than the ambient sound speed then the radio structure has 
been stable for a minimum of 3 x 10 8 yr. The surface brightness 
edge could be a weak shock generated by a merger at that time or 
more likely, given that it can be observed around most of the cluster 
core, it could be a cold front produced by the cool-core gas slosh- 
ing in the cluster potential well (Markevitch et al. 2000, Vikhlinin 
|et al.p OOl ). The swirl feature seen in the projected temperature 
map would also support this interpretation (Fig. 1 1 ). 

HI 82 1+643 shows no obvious evidence for X-ray absorb- 
ing winds from the quasar in the Chandra HETG (Fang et al.| 
2002)1, LET G JMathur et aT][2003) or in FUSE spectra ( |Oegerie| 
et al.|2 000). However, outflows could be oriented in the plane of 
the sky and therefore not evident along our line of sight. Oegerle 
|et al.| ( |2000) observed associated absorption at the redshift of 



HI 82 1+643 in the FUSE and HST spectra (see also |Bahcall et al.| 
|1992[ > and concluded that this was unlikely to be quasar intrinsic 
absorption. Instead they proposed that the broad range of ioniza- 
tion present required a multiphase absorber that is likely located 
in the central cluster galaxy. We found that the X-ray spectra for 
the central two annuli, that cover the spatial extent of the central 
galaxy, did not require any additional absorption above the Galac- 
tic column density determined by Kal berla et al.| (2005 1. However, 
these annuli contain a significant fraction of the quasar PSF and 
any additional absorption contribution to the spectrum was difficult 
to determine. We compared soft (0.5 — l.OkeV and 1.0 — 2.0keV) 
and hard band (3.0 — 7.0keV) images to see if there were any ob- 
vious potential sites for absorption. The X-ray depression at 24kpc 
from the nucleus was clearly visible in both soft band images, sur- 
rounded by cool rims of gas, however it was not visible in the hard 
band image. The hard band image contains far fewer counts and 
was dominated by the quasar PSF. Therefore a small cavity would 
be particularly difficult to distinguish. However, this could also be 
interpreted as a cool extended arm of emission, similar to those in 
the N and NW directions, with a superimposed absorption region. 

6.2 Accretion Mechanisms 

The total energy output, radiative and mechanical, of the AGN 
can be used to infer an accretion rate for H1821+643, assum- 
ing an accretion efficiency for each mechanism. Using a bolomet- 
ric correction of 50 from Vasudevan & Fabian (20071 with the 
quasar luminosity in the energy range 2— lOkeV Lx = (4.2 ± 
0.1) xlO 45 ergs 1 , we calculated the quasar bolometric luminosity 
Lt,oi ~ 2x 10 47 ergs _1 . However, this could have been underesti- 
mated if a significant number of X-ray photons from the quasar 
were scattered by the ICM. By extrapolating the electron density 
profile in to ~ 0. 1 cm~ 3 , we estimated an optical depth for the cen- 
tral 30kpc region of T ~ 0.002. Therefore, only ~ 1 in 1000 pho- 
tons from the quasar would be scattered by the ICM and this effect 
is insignificant. 

For a radiative efficiency e = 0.1, the mass accretion rate 
required to power the quasar is therefore M acc = L\, \/ec 2 ~ 
4OM yr _1 . The mechanical power calculated by cavity heating 
is several orders of magnitude lower than the quasar luminosity 
and therefore is not included in an estimate of the required accre- 
tion rate. The calculated mass accretion rate for H1821+643 can be 
compared with the theoretical Eddington and Bondi accretion rates. 

The Eddington accretion rate can be expressed as 

M E dd = 6£( m bh \ (3) 
Meyr- 1 e \3x 1O 9 M 7 

Using the Ha, Hj3 and Mgll line widths and the optical/UV 
continuu m luminosity (|Kol man et aL| |1993| >, together with the 
recipes in McG ill et al.| j2008) (their section 4.1 and references 
therein) we estimated a black hole mass M BH ~3x 10 9 M Q for 
HI 82 1+643. The uncertainties in the bolometric correction and the 
black hole mass were at least 50% and therefore subsequent cal- 
culations were assumed to be only an estimate. Using the esti- 
mate of black hole mass, the Eddington luminosity was calculated 
to be Z-Edd ~4x 10 47 ergs~' , giving an Eddington accretion rate 
M^dd ~ 70 M yr~' and a high Eddington ratio Z-bol/^Edd ~ 0.5. 

Assuming spherical symmetry and negligible angular momen- 
tum, the Bondi rate is the rate of accretion for a black hole with an 
accreting atmosphere of temperature, T, and density, n e , l [Bondi| 
|1952[ > and can be expressed as 
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^Bondi 

M Q yr _1 



-3/2 



M B h 
3xl0 9 M H 



(4) 



for an adiabatic index y = 5/3. This accretion occurs within the 
Bondi radius, rg on( ji, within which the gravitational potential of the 
black hole dominates over the thermal energy of the surrounding 

gas, 



r Bondi 

kpc 



Mm 



M B h 

3 X 10 9 Mrr, 



(5) 



The Bondi accretion rate is therefore an estimate of the rate 
of accretion from the hot ICM directly onto the black hole. The 
Bondi radius is not resolved for H1821+643 so the electron den- 
sity and temperature from the innermost annulus were used as 
lower and upper limits, respectively, on the accretion atmosphere 
properties producing an underestimate of the Bondi accretion rate. 
For an innermost temperature 2.4±0.2keV and electron density 
0.046 ± 0.002 cm -3 , the Bondi accretion rate is estimated to be 
^Bondi ~ 0.001 M yr _I at the Bondi radius r Bon( ji ~ 40pc. Al- 
though this estimate of the Bondi accretion rate is expected to lie 
below the true value, the implied accretion rate for the quasar ex- 
ceeds this by a factor of ~ 10 4 . 

Following | Allen et al.| ((2006), we also estimated the Bondi 
accretion rate by extrapolating the temperature and density of the 
cluster gas in to the Bondi radius. This was achieved by using 
powerlaw fits to the inner annuli of the DSDEPROJ deprojected 
temperature and density profiles in Fig. [T2] The steep temper- 
ature (T(r) o= r ' 6 ) and density gradients (n e (r) o= r~ 12 ) in the 
cool-core produced an upper limit on the Bondi accretion rate of 
A^Bondi ~ 6M yr~' at rgondi ~ 0.5 kpc. However, an extrapolation 
over two orders of magnitude in radius is unlikely to provide a good 
estimate of the Bondi accretion rate. Allen et al. (2006) analysed 
elliptical galaxies, which had shallow central temperature gradi- 
ents and X-ray data permitting a measurement of the gas properties 
within one order of magnitude of the Bondi radius. 

Compton cooling of the accretion flow by radiation from the 
quasar could potentially increase the flow of material from a hot 
atmosphere onto the black hole (Fabian & Crawford 1990). Since 
Compton scattering between strong UV photons and hot electrons 
can dominate the cooling process for the ICM in the cluster core, 
accretion by Compton cooling can supply gas to the central black 
hole and, as shown by Fabian & Crawford ( 1990 ), the accretion rate 
by Compton cooling flow can be larger than that of Bondi accretion 
alone. Using the cluster gas properties in the innermost annulus, 
we determined that Compton cooling will dominate bremsstrahlung 
cooling within a radius ~ 5 kpc and rapidly cool the gas down 
to the Compton temperature. 

The Compton temperature was calculated using the quasar 
spectral energy distribution (SED) from the optical to the soft X- 
ray presented in Kol man etaT| ( |1993| >. We fitted the SED in XSPEC 
with their Kerr disk model for the UV bump (Fig. 16(c) in Kol- 
|man et al.| 19 93 ) and a broken power-law to provide an extension 
up to 100 keV, as described in Vasudevan & Fabian (2009). There- 
fore, with the bulk of the emitted quasar luminosity in the UV, the 
Compton temperature is Tq ~ 0.4keV. 

If the accreting material is quickly Compton cooled within 
~ 5 kpc, as it approaches the Bondi radius, then the Bondi radius 
will grow from rB on di ~ 40pc to ~ 250pc and the Bondi accre- 
tion rate increases from MBondi ~ 0.001 Mq yr~' to ~ 6M a yr~ 1 . 



Therefore, Bondi accretion of rapidly Compton cooled material 
could provide a significant source of fuel for the quasar. 

Cold accretion mechanisms, such as the cold feedback sce- 
nario suggested by Pizzolato & Soker (2005 ), could provide an al- 
ternative route for accreting material. Within the cooling region of 
HI 82 1+643, up to 300 ± lOOM yr _1 could be cooling down be- 
low X-ray temperatures. Pizzolato & Soker ( 2005 ) proposed that 
overdense blobs of cool gas from within the entire cooling region 
rapidly cool and sink towards the central black hole or condense 
to form cold molecular clouds and stars. A large influx of material 
onto the black hole, caused by an accretion disk instability for ex- 
ample, would then trigger an outburst. For cold accretion to provide 
a source of fuel for the quasar, a significant fraction of the gas cool- 
ing out of the X-ray must be funnelled down onto the black hole. 
Such an efficient, steady accretion of material from within the cool- 
ing region seems unlikely, however the cooling material could have 
built up in the cluster centre to form a large reservoir of cool gas 
which feeds the nucleus after a triggering event. 

Alternatively, a merger scenario could have initiated the 
quasar outburst by providing a rapid influx of fuel for the central en- 
gine. The precession in the jet axes suggested by the observations of 
Blundell & Rawlings ( 2001 1 could have arisen if the central engine 
was composed of two orbiting black holes that had not yet merged. 
A merger scenario could therefore potentially have triggered the 
quasar outburst and the jet precession. However, the galaxy clus- 
ter's deep gravitational potential will tidally strip a merging sub- 
cluster of its gas in the cluster outskirts, making it difficult to chan- 
nel the merging material effectively down onto the central black 
hole. 

6.3 Impact of the quasar on the galaxy cluster core 

Fig. [T3] compares the radiative gas cooling time profiles of 
HI 82 1+643 with a sample of 16 cool-core and non-cool-core clus- 
ters with Chandra archival observations. The clusters were roughly 
selected on redshift (z = 0.06 — 0.25) and the number of X-ray 
counts in the Chandra observation. The first criterion ensured a 
similar intrinsic spatial extent on the ACIS chip to HI 82 1+643 
and the second provided an equivalent resolution with at least eight 
deprojected radial bins (~ 40,000 counts). In addition, three non- 
cool-core clusters, Abell 1650, Abell 2142 and Abell 2244, were 
included to serve as a contrast. 

HI 82 1+643 has a cooling time profile typical of a strong cool- 
core cluster with a more passive central engine, such as Abell 1835. 
Underneath the quasar, the large-scale properties of the galaxy clus- 
ter have not been strongly affected and the ICM is continuing to 
evolve in a similar way to other cool-core clusters. 

This period of quasar activity in HI 82 1+643 is unlikely to be 
long-lived. For such a high accretion rate of 40M a yr _1 , the e- 
folding time of the SMBH growth, assuming a radiative efficiency 
e = 0.1, is ~ 10 8 yr and in this time the SMBH would need to ac- 
crete ~4x 10 9 M of material. A period of rapid black hole growth 
longer than ~ 10 8 yr would therefore result in a black hole mass far 
larger than the typical upper limit observed (eg. Marco ni et al. 12004] 
and references therein). On this timescale, we would not expect ra- 
diative heating from the quasar to significantly alter the gas cooling 
time calculated in large radial bins. However, to quantify the effect 
of the active nucleus on the cluster gas we have run Cloudy (ver- 
sion 08.00) (Ferland et al. 1998 ) models appropriate to conditions 
in the regions 15 — 30 kpc (projected) and 22 — 39 kpc (DSDEPROJ 
deprojected) from the nucleus. 

The Cloudy simulations are described in detail in Appendix B. 
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Figure 15. DSDEPROJ deprojected gas cooling times for a sample of cool- 
core and non-cool-core galaxy clusters. 



In summary, we generated the intrinsic continuum shape for the ac- 
tive nucleus using the quasar SED model described in section |6T2] 
A model of the cluster gas, with suitable temperature 2.4keV and 
metallicity 0.4 Zq, at the distance of the extracted X-ray spectrum 
was then irradiated with this quasar continuum. The emitted spec- 
trum of the Cloudy model, together with a model accounting for 
the quasar PSF, were fitted to the observed spectra from the pro- 
jected (15 — 30 kpc) and deprojected (22 — 39kpc) regions. We also 
produced a Cloudy model for only the cluster emission to provide 
a simple comparison for the Cloudy model with the added quasar 
contribution. 

The addition of the Cloudy quasar component did not signif- 
icantly alter the shape of the model spectrum (Fig. |16|l. However 
in this model, photoionization of the cluster gas produced a signif- 
icantly improved fit to the Fe K a line blend at 5 keV in the pro- 
jected spectrum (x 2 = 130 reduced to x 2 = 1 10 for 100 degrees of 
freedom). This improvement might be due to continuum pumping 
( |Porter & Ferland|2007| >. 

This could provide an explanation for the sudden drop in 
metallicity seen between 15 and 26 kpc in the projected profile (Fig. 
1 1 0| > . This effect will obviously be stronger closer to the nucleus, 
therefore it is unlikely to be observed in the deprojected metallicity 
profile which did not extend as close in to the quasar as the pro- 
jected profile. 

However, whilst our Cloudy simulations suggest that the 
quasar could be significantly photoionizing the gas close to the nu- 
cleus, there are several other possible explanations for the observed 
drop in core metallicity. Buote & Canizares (1994) argue that fit- 
ting intrinsically multitemperature spectra with single-temperature 
models artificially produces a lower metallicity value. There could 
be a significant contribution from a second temperature component 
in the observed projected spectrum of H1821+643, such as a hot- 
ter projected cluster emission or emission from cooling gas in the 
cluster core. 

We therefore added a second tempera- 
ture component to our XSPEC model to give a 
PHABS(POWERLAW)+PHABS(MEKAL+MEKAL) model which 
was fitted to the projected spectrum (15 — 30kpc). The second 



Figure 16. The cluster projected spectrum from 15 — 30kpc is shown with 
the best-fit Cloudy and XSPEC models superimposed. Upper: the Cloudy 
models for the cluster gas (black solid line) and the cluster gas irradiated 
by the quasar (red dashed line). Lower: the Cloudy model irradiated by 
the quasar (red dashed line) and the two-temperature MEKAL model with a 
5.0keV (blue dot-dashed line) or a 0.5keV (purple solid line) second tem- 
perature component. The contribution of the quasar PSF was included as a 
PHABS(POWERLAW) component in each model. The three models produce 
a very similar spectral fit. 



temperature component was first fixed to 5.0keV and the metal- 
licity in both temperature components was fixed to 0.5Z e . Fig. 
[Trj| (lower panel) shows that the addition of a second temperature 
component produced an improved fit to the Fe K a line blend. 
However, a cooler second temperature component of 0.5 keV was 
found to provide a significantly better fit to the projected spectrum 
than the 5.0keV component (x 2 = 94 for 0.5 keV compared 



to x = 109 for 5.0keV with 99 degrees of freedom). Fig. 16 
suggests that a second temperature component, hotter or colder, or 
photoionization of the cluster gas could explain the poor fit to the 
Fe K a line blend. 

Assuming pressure equilibrium, the volume filling fraction of 
the 0.5 keV temperature component i is 



fi- 



£jT 2 



(6) 



where £,- is the emission measure (XSPEC normalization) of the 
component i, T, is the temperature and j sums over all the com- 
ponents. The 0.5 keV temperature component has a volume fill- 
ing fraction of 0.2%, which equates to a gas mass of M gas ~ 
3 x 1O 9 M . With a very short radiative gas cooling time of only 
50Myr, the 0.5 keV temperature component could produce a mass 
deposition rate of ~ 5OM yr~ 1 in the central 30kpc in the absence 
of heating. 

Finally, an additional power-law component from inverse- 
Compton emission could also provide an explanation for the poor 
fit to FeKa. A power-law component would raise the level of 
the continuum emission in the projected spectrum and reduce the 
prominence of the emission lines, in a similar way to the extra 5 keV 
temperature component. However, coincident cluster emission and 
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quasar PSF made it difficult to put any constraints on inverse- 
Compton emission and therefore we cannot determine whether this 
would resolve the model fit. 

Although from our results it is not possible to determine a pre- 
ferred explanation for the poor fit to the Fe K a line blend, we ar- 
gue that a combination of the above would be most likely. The pro- 
jected spectrum will contain a contribution of hotter (5.0— 8.0keV) 
gas from the outer cluster layers. Our Cloudy simulations suggest 
that the quasar emission should produce a significant amount of 
photoionization in the central 30kpc. Observations of strong cool- 
core clusters with the XMM-Newton Reflection Grating Spectrom- 
eter have confirmed that many have a second cooler temperature 
component in the central region (Peterson et al. 2003 ; Kaastra et al. 
[2004l|Sanders et al.|2008| . Finally, although we are unable to quan- 
tify any inverse-Compton emission, this could also be a factor close 
in to the nucleus. 



6.4 AGN Unified Accretion Model 



Churazov etaEl <|2"005 ) (see also Di Matteo et al. 2005 ; Croto n et al] 



2006| > argue that a parallel can be drawn between the evolution of 
Galactic blac k holes (eg. |Fender et al.||1999| |Gallo et al.||2003^ 
and SMBHs ( [Merloni et al.||2003||Maccarone et al.|2003[ |Falcke| 
|et al.|2004[ [Chiaberge et al. 2005}. SMBHs and their surrounding 
medium are therefore expected to evolve through two stages. Early 
on, the SMBH grows rapidly by accreting cooling gas near the Ed- 
dington rate with weak feedback on the surrounding medium but 
a very high radiation efficiency (eg. |Yu &"Tr emaine 2002). This is 
known as the 'quasar-mode' which terminates when, despite the 
weak coupling to the gas, the SMBH produces enough heat to 
suppress cooling. With a lower accretion rate, the system moves 
to a more passively evolving stage, 'radio-mode', where efficient 
mechanical feedback sustains the ICM at X-ray temperatures. Re- 
cent observations of nearby low-luminosity SMBHs in clusters and 
galaxies suggest that mechanical feedback from the AGN, in the 
form of bubbles, shocks and sound waves, is effective at heating 
the cluster g as (eg.|Fabian et al.|2003||Birzan et al.|2004||Dunn &| 



|Fabian! 2006, Raf ferty et aT[2 006; Fabi an et al.|2006 



With rapid accretion at half of the Eddington rate, HI 82 1+643 
is currently in 'quasar-mode', corresponding to the high/soft state 
of black hole X-ray binaries where a geometrically thin, opti- 
cally thick accretion disk forms (Shak ura & Sunyaev|1973} . How- 
ever, quantifying the quasar feedback into the surrounding ICM in 
HI 82 1+643 is difficult. HI 82 1+643 has a gas cooling time profile 
typical of a strong cool-core cluster with a passive central engine 
indicating that the quasar has not strongly affected the large-scale 
cluster gas properties. The Cloudy simulations detailed in section 
|6.3| also suggest that the radiative output of the luminous quasar has 
little impact on the bulk of the cluster gas. There is some indication 
of mechanical heating, although this has a low heating power and 
would only currently compensate for ~ 10% of the ICM cooling 
losses. As discussed in section |6~T) there is no obvious evidence 
for X-ray absorbing winds in HI 82 1+643 in the existing observa- 
tions. However, outflows might not be evident along our line of 
sight and therefore cannot be ruled out. Although the cluster gas 
properties do not indicate any efficient feedback from quasar winds 
on the scales considered (r > 15kpc) this does not exclude an effi- 
cient coupling to the cold galaxy gas (see Chartas et al. 2002, 2003 
[20071 |Pounds et al.|2003||Reeves et al.|2009| >. 

HI 82 1+643 could provide an example of the early stage of 
evolution for these systems where the SMBH has not evolved be- 
yond 'quasar-mode' and cannot produce enough heat to suppress 



cooling in this rich cluster. The underlying cluster may exhibit the 
early cool-core environment before 'radio-mode' heating switches 
on to efficiently reheat the gas. It seems implausible that the quasar 
could remain switched on for longer than ~ 10 8 yr because in that 
time the black hole mass will increase to ~ 1O 1O M0, greater than 
the typical observed upper limit (eg. Marcon i et al.|2004| >. 

However, the FR I radio structure shown by HI 82 1+643 ap- 
pears to be at odds with the simple picture of 'quasar-mode' feed- 
back. Around 10% of quasars are known to have strong radio emis- 
sion (eg. |Ivezic et al.|2 002 l indicating that SMBHs with high ac- 
cretion rates can be associated with powerful outflows. By draw- 
ing comparisons with observations of microquasars, Nipot i et al.| 
(2005) suggest a parallel between radio-loudness in quasars and 
the flaring mode in microquasars, where radio-loudness is simply 
defined as extended radio emission (FR I or FR II). Microquasars 
have been observed to have two distinct modes of energy output: 
one producing roughly steady radio emission that is coupled with 
X-ray emission (coupled mode), the other producing strongly vari- 
able, flaring radio emission that appears to be decoupled from the 
X-ray emission (flaring mode). Nipo ti et al.| ( [2005] > further asso- 
ciate radio-quiet quasars with non-flaring/coupled states of micro- 
quasars, in either low/hard or high/soft states. 

In this picture, H1821+643 would be classed as a radio-loud 
quasar in flaring mode, where the nucleus is producing the extended 
FR I radio structure. These flares could be analogous to the short, 
radio outbursts that have been observed in transient X-ray binaries 
during the transition from the low/hard to high/soft state when the 
accretion rate approaches the Eddington limit ( |Fender et al.|2004| >. 
However, whilst black hole accretion can be scaled from micro- 
quasars and X-ray binaries up to quasars, the environments of these 
objects vary significantly. It is unclear how the rich cluster environ- 
ments and accretion history of quasars will affect the conclusions 
drawn from smaller scale systems. 



7 CONCLUSIONS 

By accurately subtracting the quasar contribution, we were able 
to determine the properties of the surrounding ICM down to 3" 
from the quasar. The temperature of the cluster gas decreases 
from 9.0±0.5keV beyond ~ 200 kpc down to 1.3±0.2keV inside 
20kpc, with a short central radiative cooling time of 1 .0 ± 0. 1 Gyr, 
typical of a strong cool-core cluster. By comparing the cooling time 
with a sample of cool-core clusters with 'radio-mode' AGN, we de- 
termined that the quasar does not appear to have had a significant 
impact on the large-scale cluster gas properties. 

The cluster core features several extended arms of emission, 
which could have been dragged out by the expanding FR I struc- 
ture. However, we found that the radio emission is not clearly re- 
lated to the X-ray gas morphology. The nature of the surface bright- 
ness edge at ~ 15" could not be unequivocably determined from 
the temperature and density profiles. This could be interpreted as 
a weak shock generated by the expansion of the radio lobes, or a 
cold front produced by the cool-core sloshing in the cluster's grav- 
itational potential. There is some evidence for cavity heating in the 
cluster core which could currently compensate for ~ 10% of the 
ICM cooling losses. 

The derived accretion rate of 4OM yr _1 is approximately 
half of the quasar's Eddington limit. We showed that Bondi accre- 
tion, boosted by Compton cooling of the accretion material by the 
quasar radiation, could provide a significant proportion of the re- 
quired fuel. However, without resolving the cluster gas properties 
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at the Bondi radius, we could not determine whether this mecha- 
nism could supply all of the required 4OM yr and potentially 
provide a self-sustaining fuel source for the quasar. We have also 
discussed the alternative mechanisms of cold accretion or a rapid 
infall of a large quantity cool gas provided by a subcluster merger. 

In the context of AGN unified accretion models, we suggest 
that HI 82 1+643 could provide an example of the early stage of evo- 
lution for these systems where the SMBH has not evolved beyond 
'quasar-mode' and cannot produce enough heat to suppress cool- 
ing. The surrounding rich cluster could therefore exhibit the early 
cool-core environment before the onset of 'radio-mode' heating ef- 
ficiently reheats the gas. However, the powerful FR I radio struc- 
ture observed in H1821+643, and in around 10% of quasars (eg. 
|Ivezic et al.|2002| >, conflicts with this simple scenario of 'quasar- 
mode' feedback. Following the comparisons with observations of 
microquasars suggested by Nipoti et al. ( 2005 I, we consider that 
H1821+643 is in a flaring mode where the radio emission is de- 
coupled from the X-ray emission. These flares could be analogous 
to the short, radio bursts observed in transient X-ray binaries dur- 
ing the transition from the low/hard to high/soft state (Fender et al. 
|2004| >. However, it is unclear how the conclusions drawn from the 
smaller scale systems will be affected by the rich cluster environ- 
ment and the quasar accretion history. 
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APPENDIX A: SIMULATING THE QUASAR PSF 
Al Readout streak spectrum 

For heavily piled up sources, M. Bautz suggested using the ACIS 
readout streak to extract a source spectrum (Gaetz 2004). In this 
observation, ACIS accumulates events for a frame exposure time of 
3.1s and then reads out the frame at a parallel transfer rate of 40/1 s , 
taking a total of 0.04104 s to read out the entire frame. During read- 
out, the CCD still accumulates events, which are distributed along 
the whole column as it is read out, creating a continuous streak for 
very bright point sources. 

We extracted a spectrum in the energy range 0.5 — 7.0keV 
from two narrow regions (300 x 8 pixel 2 ) of the readout streak po- 
sitioned either side of the central source and avoiding as much of 
the cluster emission as possible. Any remaining cluster emission in 
the spectrum was then subtracted by using background regions ad- 
jacent to the readout streak. Using a wider extraction region for the 
readout streak was found to increase the uncertainty in the spectrum 
without altering the results. 

The accumulated exposure time in the streak during frame 
transfer was (300 x 2) x 27463 x 40/ts = 659.115, where 27463 
was the total number of frames. Therefore, the readout spectrum 
count rate was multiplied by a factor of 129.17, the ratio of the ob- 
servation exposure time to the readout streak exposure time. Pre- 
vious analyses of Chandra readout streak data have indicated that 
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Figure 17. HI 82 1+643 readout streak spectrum and spectral fit. The spec- 
tram has been scaled by the ratio of the observation exposure time to the 
effective streak exposure time. 



a gain correction of +2.5% for the ACIS-S3 chip (Gaetz 2004). 
Frame transfer is effectively a different ACIS mode for which there 
is likely to be a different response to arriving events compared to 
the calibrate timed exposure mode. However, for the HI 82 1+643 
readout streak on the ACIS-S3 chip, we found that the recom- 
mended +2.5% gain correction did not significantly alter the spec- 
tral fitting results and have not applied it to our analysis. 

The response files for the transfer streak spectrum were gen- 
erated using the standard CIAO tools MKACISRMF and MKWARF. 
The response matrix (RMF) was extracted over the readout streak 
region used in this analysis. However, the effective area (ARF) was 
determined using a region at the location of the quasar image. The 
X-rays in the transfer streak actually hit the detector at the direct 
image, and their detected position is an artifact of the chip readout. 
The scaled readout streak spectrum was grouped with a minimum 
of 50 counts per spectral bin, corrected for the applied count rate 
scaling factor. 

The scaled readout streak spectrum, shown in Fig. [17] was 
fitted with an absorbed power-law model in the spectral fitting 
package XSPEC version 12. The best fit photon index was con- 
sistent with the Fang et al. (2002jl analysis of the Chandra HEG 
and MEG spectra and the previous ASCA result ( Yamashita et al. 
|1997| >. Therefore, we fixed the photon i ndex to the more accu- 
rate HEG/MEG result, T = 1.761" 1 """ 



( [Fang et al.||2002| . The 
Galactic absorption parameter was left free, giving a value of 
0.003 x 10 22 cm~ 2 



with an upper limit of 0.04 x 10 22 cm 2 



Fang 



et al. (2002) inferred from the low best-fit column density ;ih that 
H1821+643 has a soft excess. However, the best-fit absorbed pow- 
erlaw is a sufficiently good description of the readout streak spec- 
trum so we have used this simplistic model in our analysis. 

The best-fitting model gave an unabsorbed flux of 

in the energy range 2 — lOkeV 



1.45 



-0.04 
-0.03 



x 10 11 erg cm 2 s 



and a reduced Xv = 1-6- The quasar flux measured from the 
readout streak spectrum falls between the values determined for 
this object by pang et al] ( |2002^ and |Yamashita et al.| ( fT997) . This 
gives a 2 — lOkeV quasar luminosity of (4.2 ± 0. 1 ) x 10 45 erg s~ 1 . 

This spectral model from the quasar readout streak was used 
as the principal component in the ChaRT simulation of the quasar. 
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A2 ChaRT simulations 

The Chandra Ray Tracer (ChaRT), in conjunction with the MARX 
software, was used to simulate the Chandra PSF produced by the 
quasar. Following the ChaRT analysis thread^] the readout streak 
spectrum and observation exposure time were used to generate a 
ray-tracing simulation. The output from ChaRT was then supplied 
to MARX version 4.4.0 which projected the ray traces onto the de- 
tector and applied the detector response. 

ChaRT provides a user-friendly interface to the SAOTrace 
semi-empirical model ( [Jerius et' al. 1995) which is based on the 
measured characteristics of the mirrors, support structures and baf- 
fles. The model is then calibrated by comparison with actual ob- 
servations. However, while the PSF core and inner wing region 
match well with observations ( Jerius 2002), the PSF wings be- 
yond 10", produced by mirror scattering, seem to be underpredicted 
by the raytrace model. At the moment, the SAOTrace model does 
not model the dither motion of the telescope or include residual 
blur from aspect reconstruction errors. MARX includes the effects 
of these errors using the DitherBlur parameter. We have therefore 
tested the use of SAOTrace to simulate the H1821+643 quasar spec- 
tra by performing an identical analysis on the quasar 3C 273, which 
has no surrounding cluster emission. Both of these sources are close 
to on-axis and only the PSF core and inner wing region are of inter- 
est to our analysis. Beyond 10" the cluster emission in H1821+643 
will dominate and an accurate model of the quasar PSF becomes 
much less important. 

A2.1 Test object 3C 273 

3C273 is a bright, nearby (z = 0.158) radio-loud quasar with a 
jet showing superluminal motion and has been intensively studied 
at different wavelengths (eg. |Courvoisier|1998[|Soldi et al.|2008^ . 
3C 273 was observed with Chandra ACIS-S for 160 ks to study the 
bright X-ray jet in detail ( |Jester et al.|2 006). Here we used one of 
the 40 ks Chandra ACIS-S exposures of this object (obs. id 4879) to 
test whether a ChaRT simulation using the readout streak spectrum 
would recover the correct result. This observation of 3C 273 was 
similar to that of HI 82 1+643, featuring a piled up point source, 
strong readout streak, bright PSF wings but crucially no detected 
extended emission (the region containing the jet was excluded). 
Therefore this observation of 3C 273 allowed us to test our method 
of determining the quasar spectrum in a series of annuli and com- 
pare the result with the observed spectra without the added compli- 
cation of superimposed cluster emission. 

As in section Al for H1821+643, the readout streak spectrum 
was extracted using two regions of the transfer streak either side 
of the nucleus. The spectrum was fitted with an absorbed power- 
law model which provided a reasonable fit for our purposes over 
the required energy range, reduced Xv = 1-2- This simple model 
gave a best-fitting photon index T = 1.64^q ^, Galactic absorp- 
tion upper limit of «h < 0.015 x 10 22 cm~ 2 and unabsorbed flux 
^6.5^Q 2 j x 10~ n ergcm~ 2 s~' in the 2— lOkeV energy range 
(Fig. [T8| >. These parameters fall within the expected range defined 
by the long term variability of this source (Chernyakov a et al.|2007| 
|Soldi et al.|20 08 and references therein). Although 3C 273 has also 
been found to have a soft excess extending up to 2keV {Turner etal.| 
[l985l|Staubertetal.|1992) we have ignored this for our simplistic 
analysis. 3C 273 has a slightly lower photon index than H 1 82 1 +643 
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Figure 18. 3C 273 readout streak spectrum and spectral fit. The spectrum 
has been scaled by the ratio of the observation exposure time to the effective 
streak exposure time. 



which produces a wider PSF because the PSF is energy-dependent 
(Figure]?). 

The best-fitting model to the readout streak spectrum was used 
as the main input for the ChaRT simulations of 3C273. Spectra 
were extracted in two annuli, excluding the readout streak, with an 
innermost radius of 2.5" outside of which pileup was determined 
to be negligible. The exposure time of each frame in this obser- 
vation of 3C273 was 0.4s, compared with 3.1s for H1821+643, 
which greatly reduced the number of piled up photons. These an- 
nuli were selected to contain the required minimum of 3000 counts 
for a good deprojection and to cover a similar area of the PSF wings 
to the regions used in the analysis of H1821+643. The observed and 
simulated spectra were deprojected and are shown overlaid in Fig. 
[19| The simulated spectra were generally well-matched to the ob- 
served spectra. The mismatches can be mostly accounted for by 
comparison with the model of the readout streak spectrum, partic- 
ularly at ~ 3keV where the model was a significant overestimate. 
The best-fitting parameters for an absorbed power-law model fitted 
to the observed and simulated deprojected spectra were also found 
to be consistent within the errors. 

A ChaRT simulation based on the readout streak spectrum was 
therefore expected to provide an adequate description of the quasar 
spectrum for the series of radial bins used in this analysis. 



A2.2 ChaRT Simulation of HI 821 +643 

A ChaRT simulation of the quasar HI 82 1+643 was generated for 
the same position on the detector using the model of the readout 
streak spectrum. ChaRT was run for an effective exposure time 
three times longer than the total observation time to sample the 
range of possible optical paths in the HRMA and reduce statisti- 
cal errors. The raytraces were then projected onto the detector with 
MARX to create an event file. Simulated quasar spectra were ex- 
tracted from the MARX event file and combined with appropriate 
response files for the simulatiorj^] 

The width of the simulated PSF depends on the MARX 



6 See http://cxc.harvard.edu/chart/ 
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7 See http://space.mit.edu/CXC/MARX/examples/ciao.html 
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trum is caused by an incorrect calculation of the effective area. This 
also causes the variation in these parameters between the radial 
bins. The effective area correction for each spectrum, provided by 
the ancillary response files (ARFs) assumes that all of the PSF falls 
within that extraction region. This was true for the cluster spectrum 
in each annulus but not for the quasar spectrum, where only the PSF 
wings were included (less than 5% of the total PSF). The effective 
area correction for the quasar spectrum should therefore be multi- 
plied by the encircled energy fraction, which varies as a function 
of energy. When the ARF spectral response was multiplied by the 
EEF calculated in a series of broad energy bands from the ChaRT 
simulation, the input quasar parameters are recovered. 

We confirmed our earlier analysis that pileup was not impor- 
tant beyond a radius of 3" from the quasar by analysing a MARX 
simulation including pileup. The pileup simulation produced es- 
sentially identical results to the simulations without pileup and all 
spectral parameters were consistent within errors. 
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Figure 19. Comparison of the observed deprojected spectra (black circles) 
and ChaRT simulated deprojected spectra (red triangles) for two annuli in 
3C 273. Upper: 2.5 - 4.9" Lower: 4.9 - 9.8". 

DitherBlur parameter, which is a statistical term combining the as- 
pect reconstruction error, ACIS pixelization and pipeline pixel ran- 
domization. The default MARX value is 0.35", however the magni- 
tude of the aspect blur is observation dependent. The MARX simu- 
lation of H1821+643 was repeated for DitherBlur values of 0.25", 
0.30" and 0.40" but for all three cases the quasar-subtracted cluster 
surface brightness profile was consistent within statistical errors. In 
addition, all spectral parameters were found to be consistent within 
errors to the 0.35" simulation. We therefore proceeded with the 
DitherBlur parameter set to its default. 

The simulated quasar spectra were extracted in identical re- 
gions to regions 1-11 (projected) and regions a-b (deprojected) 
shown in Fig. [5] The simulated quasar spectra from regions a-b 
were deprojected with DSDEPROJ, as detailed in section |4~3) to en- 
sure they have been processed in an identical way to the observed 
deprojected spectra. By determining the quasar spectral parameters 
from the simulation for the same regions, we can fix these values in 
the combined fit to the observation of the quasar and cluster. 

The simulated quasar spectrum from each region was fit- 
ted with an absorbed power-law model PHABS(POWERLAW) in 
XSPEC; the best-fitting parameters are given in table[T| 

The difference between the best-fit n# and power-law param- 
eters (Table[TJ and their input values from the readout streak spec- 



APPENDIX B: CLOUDY SIMULATIONS 

The most obvious way to model the intracluster gas using Cloudy 
is with the CORONAL EQUILIBRIUM command. However, this sets 
up a constant temperature model which is not appropriate when we 
want to consider the effect of the ionizing flux from the quasar on 
this gas. Also, it is not clear that a pure cooling model is actually 
physically appropriate since the cores of clusters are known to be 
deficient in cool gas and this is widely thought to be due to the 
presence of a heat source (eg. Peterson & Fabian 2006; McNamara 
|& Nulsen|2007| >. We therefore chose to model the cluster emission 
by using an arbitrary heat source and adjusting the heating rate to 
obtain the observed temperature of the cluster through a balance of 
heating and cooling. This was achieved using the HEXTRA com- 
mand. 

Therefore, we first simulated cluster emission by creating a 
single zone model heated using only the HEXTRA pure heating 
command. Metal abundances were set to match those measured 
by fitting a single temperature mekal model to the X-ray spec- 
tra, 0.4 Z . The heating rate was adjusted until the temperature 
of the gas matched that measured in the cluster in these regions 
(15 - 30kpc and 22- 39kpc), 2.78 x 10 7 K (2.4 keV). This model 
spectrum was then renormalized to fit the cluster spectrum and 
the chi-square statistic evaluated. This Cloudy model for the clus- 
ter emission provided a simple comparison for the Cloudy model 
with the added quasar contribution; any differences could then be 
attributed to the quasar input and not a discrepancy between the 
Cloudy and MEKAL models of cluster emission. 

To add in the quasar contribution we first generated the in- 
trinsic continuum shape for the active nucleus using the quasar 
SED model described in section [672] To approximate the correct 
spectrum to irradiate the cluster gas at the distance of our ex- 
tracted X-ray spectra we set up a Cloudy model with this contin- 
uum and allowed it to pass though cluster gas at a temperature of 
2.78 x 10 7 K (2.4keV) for an equivalent hydrogen column density 
of projected 2.3 x 10 21 cm~ 2 (15kpc x 0.05 cm -3 ) or deprojected 
3.4 x 10 21 cirr 2 (22kpcx0.05cnr 3 ). 

The transmitted spectrum from the output of this model was 
then added to the pure cluster simulation, allowing the HEXTRA 
component to vary so that the observed temperature of the cluster 
was still reproduced, 2.78 x 10 7 K (2.4keV). The emitted spectrum 
from this Cloudy model together with a model accounting for the 
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quasar PSF (section |4~T| l were fitted to the observed spectra from 
the projected (15 — 30kpc) and deprojected (22 — 39kpc) regions. 



REFERENCES 

Allen S. W., Dunn R. J. H., Fabian A. C, Taylor G. B., Reynolds 

C. S., 2006, MNRAS, 372, 21 

Allen S. W., Ettori S., Fabian A. C, 2001, MNRAS, 324, 877 
Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 53, 
197 

Arnaud K. A., 1996, in Astronomical Society of the Pacific Con- 
ference Series, Vol. 101, Jacoby G. H., Barnes J., ed, Astronom- 
ical Data Analysis Software and Systems V, p. 17 
Bahcall J. N., Jannuzi B. T., Schneider D. P., Hartig G. F, Green 

R. F, 1992, ApJ, 397, 68 
Balucinska-Church M., McCammon D., 1992, ApJ, 400, 699 
Birzan L., Rafferty D. A., McNamara B. R., Wise M. W., Nulsen 

P. E. J., 2004, ApJ, 607, 800 
Blundell K. M., Rawlings S., 2001, ApJ, 562, L5 
Boehringer H., Voges W., Fabian A. C, Edge A. C, Neumann 

D. M., 1993, MNRAS, 264, L25 
Bondi H., 1952, MNRAS, 112, 195 

Buote D. A., Canizares C. R., 1994, ApJ, 427, 86 
Carter C, Karovska M., Jerius D., Glotfelty K., Beikman S., 2003, 
in Astronomical Society of the Pacific Conference Series, Vol. 
295, Payne H. E., Jedrzejewski R. I., Hook R. N, ed, Astronom- 
ical Data Analysis Software and Systems XII, p. 477 
Chartas G., Brandt W. N., Gallagher S. C, 2003, ApJ, 595, 85 
Chartas G., Brandt W. N., Gallagher S. C, Garmire G. P., 2002, 
ApJ, 579, 169 

Chartas G., Brandt W. N, Gallagher S. C, Proga D., 2007, AJ, 
133, 1849 

Chernyakova M. et al., 2007, A&A, 465, 147 

Chiaberge M., Capetti A., Macchetto F. D., 2005, ApJ, 625, 716 

Churazov E., Sazonov S., Sunyaev R., Forman W., Jones C, 

Bohringer H., 2005, MNRAS, 363, L91 
Clarke T. E., Blanton E. L., Sarazin C. L., 2004, ApJ, 616, 178 
Courvoisier T. J.-L., 1998, A&A Rev., 9, 1 
Croton D. J. et al„ 2006, MNRAS, 365, 1 1 
Davis J. E., 2001, ApJ, 562, 575 

Di Matteo T., Springel V, Hernquist L., 2005, Nat, 433, 604 
Dunn R. J. H., Fabian A. C, 2004, MNRAS, 355, 862 
Dunn R. J. H., Fabian A. C, 2006, MNRAS, 373, 959 
Edge A. C, Frayer D. T., 2003, ApJ, 594, L13 
Fabian A. C, 1994, ARA&A, 32, 277 
Fabian A. C, Crawford C. S., 1990, MNRAS, 247, 439 
Fabian A. C, Hu E. M., Cowie L. L., Grindlay J., 1981, ApJ, 248, 
47 

Fabian A. C, Sanders J. S., Allen S. W., Crawford C. S., Iwa- 

sawa K., Johnstone R. M., Schmidt R. W., Taylor G. B., 2003, 

MNRAS, 344, L43 
Fabian A. C. et al., 2000, MNRAS, 318, L65 
Fabian A. C, Sanders J. S., Taylor G. B., Allen S. W., Crawford 

C. S., Johnstone R. M., Iwasawa K., 2006, MNRAS, 366, 417 
Falcke H., Kording E., Markoff S., 2004, A&A, 414, 895 
Fang T, Davis D. S., Lee J. C, Marshall H. L., Bryan G. L., 

Canizares C. R., 2002, ApJ, 565, 86 
Fender R. et al., 1999, ApJ, 519, LI 65 

Fender R. P., Belloni T. M., Gallo E., 2004, MNRAS, 355, 1 105 
Ferland G. J., Korista K. T, Verner D. A., Ferguson J. W., King- 
don J. B., Verner E. M., 1998, PASP, 1 10, 761 



Forman W. et al., 2005, ApJ, 635, 894 
Fried J. W., 1998, A&A, 331, L73 

Gaetz T. J., 2004, Analysis of the Chandra On-Orbit PSF Wings, 

http://cxc.harvard.edu/cal/Hrma/psf/ 
Gallo E., Fender R. P., Pooley G. G., 2003, MNRAS, 344, 60 
Hall P. B., Ellingson E„ Green R. F, 1997, AJ, 1 13, 1 179 
Hatch N. A., Crawford C. S., Fabian A. C, 2007, MNRAS, 380, 

33 

Hicks A. K., Mushotzky R., 2005, ApJ, 635, L9 
Hutchings J. B., Neff S. G, 1991, AJ, 101, 2001 
Ivezic Z. et al., 2002, AJ, 124, 2364 

Jerius D., 2002, Comparison of on-axis Chandra 
Observations of AR Lac to SAOsac Simulations, 
http://cxc.harvard.edu/cal/Hrma/psf/ARLac-onaxis.ps 

Jerius D., Freeman M., Gaetz T, Hughes J. P., Podgorski W., 
1995, in Astronomical Society of the Pacific Conference Series, 
Vol. 77, Shaw R. A., Payne H. E., Hayes J. J. E., ed, Astronomi- 
cal Data Analysis Software and Systems IV, p. 357 

Jester S., Harris D. E., Marshall H. L., Meisenheimer K., 2006, 
ApJ, 648, 900 

Johnstone R. M., Fabian A. C, Nulsen P. E. J., 1987, MNRAS, 
224, 75 

Jones T. W., De Young D. S., 2005, ApJ, 624, 586 
Kaastra J. S., 1992, in Internal SRON-Leiden Report, updated ver- 
sion 2.0 

Kaastra J. S. et al„ 2004, A&A, 413, 415 

Kalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Ba- 
jaja E„ Morras R., Poppel W. G. L., 2005, A&A, 440, 775 
Kolman M., Halpera J. P., Shrader C. R., Filippenko A. V, Fink 

H. H., Schaeidt S. G., 1993, ApJ, 402, 514 
Kriss G. A., Cioffi D. F, Canizares C. R., 1983, ApJ, 272, 439 
Lacy M., Rawlings S., Hill G. J., 1992, MNRAS, 258, 828 
Liedahl D. A., Osterheld A. L., Goldstein W. H, 1995, ApJ, 438, 
L115 

Maccarone T. J., Gallo E., Fender R., 2003, MNRAS, 345, L19 
Marconi A., Risaliti G., Gilli R., Hunt L. K, Maiolino R., Salvati 

M., 2004, MNRAS, 351, 169 
Markevitch M. et al, 2000, ApJ, 541, 542 
Markevitch M., Vikhlinin A., 2007, Phys. Rep., 443, 1 
Mathur S., Weinberg D. H., Chen X., 2003, ApJ, 582, 82 
McGill K. L„ Woo J.-H., Treu T, Malkan M. A., 2008, ApJ, 673, 

703 

McNamara B. R., Nulsen P. E. J., 2007, ARA&A, 45, 1 17 
McNamara B. R. et al., 2006, ApJ, 648, 164 
McNamara B. R. et al., 2000, ApJ, 534, L135 
McNamara B. R. et al, 2001, ApJ, 562, L149 
Merloni A., Heinz S., di Matteo T, 2003, MNRAS, 345, 1057 
Mewe R., Gronenschild E. H. B. M., van den Oord G. H. J., 1985, 
A&AS, 62, 197 

Mewe R., Lemen J. R., van den Oord G. H. J., 1986, A&AS, 65, 
511 

Navarro J. F, Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 

Nipoti C, Blundell K. M., Binney J., 2005, MNRAS, 361, 633 

O'Dea C. P. et al., 2008, ApJ, 681, 1035 

Oegerle W. R. et al, 2000, ApJ, 538, L23 

Peterson J. R., Fabian A. C, 2006, Phys. Rep., 427, 1 

Peterson J. R., Kahn S. M., Paerels F. B. S., Kaastra J. S., Tamura 

T, Bleeker J. A. M., Ferrigno C, Jernigan J. G., 2003, ApJ, 590, 

207 

Pizzolato F, Soker N, 2005, ApJ, 632, 821 
Porter R. L., Ferland G. J., 2007, ApJ, 664, 586 



© 0000 RAS, MNRAS 000, 000-000 



20 H.R. Russell et al. 



Pounds K. A., King A. R., Page K. L., O'Brien P. T., 2003, MN- 
RAS, 346, 1025 

Rafferty D. A., McNamara B. R., Nulsen P. E. J., Wise M. W., 

2006, ApJ, 652, 216 
Reeves J. N. et al., 2009, ApJ, 701, 493 

Russell H. R., Sanders J. S., Fabian A. C, 2008, MNRAS, 390, 
1207 

Sanders J. S., 2006, MNRAS, 371, 829 

Sanders J. S., Fabian A. C, 2007, MNRAS, 381, 1381 

Sanders J. S., Fabian A. C, Allen S. W., Morris R. G., Graham J., 

Johnstone R. M., 2008, MNRAS, 385, 1186 
Sanders J. S., Fabian A. C, Taylor G. B., 2005, MNRAS, 356, 

1022 

Saxton R. D., Barstow M. A., Turner M. J. L., Williams O. R., 

Stewart G. C, Kii T, 1997, MNRAS, 289, 196 
Schindler S., Castillo-Morales A., De Filippis E., Schwope A., 

Wambsganss J., 2001, A&A, 376, L27 
Schmidt R. W., Allen S. W., Fabian A. C, 2001, MNRAS, 327, 

1057 

Schneider D. P., Bahcall J. N., Gunn J. E., Dressier A., 1992, AJ, 
103, 1047 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Soldi S. et al., 2008, A&A, 486, 411 

Staubert R., Fink H., Courvoisier T. J.-L., Ulrich M.-H., Brun- 
ner H., Zimmermann U., Kendziorra E., Otterbein K., 1992, in 
American Institute of Physics Conference Series, Vol. 254, Holt 
S. S., Neff S. G., Urry C. M., ed, American Institute of Physics 
Conference Series, p. 366 
Turner M. J. L., Courvoisier T, Staubert R., Molteni D., Trumper 

J., 1985, Space Science Reviews, 40, 623 
Vasudevan R. V., Fabian A. C, 2007, MNRAS, 381, 1235 
Vasudevan R. V., Fabian A. C, 2009, MNRAS, 392, 1 124 
Vikhlinin A., Markevitch M., Murray S. S., 2001, ApJ, 551, 160 
Voigt L. M., Fabian A. C, 2006, MNRAS, 368, 518 
White D. A., Jones C, Forman W., 1997, MNRAS, 292, 419 
Yamashita A., Matsumoto C, Ishida M., Inoue H., Kii T, Mak- 

ishima K., Takahashi T, Tashiro M., 1997, ApJ, 486, 763 
Yu Q., Tremaine S., 2002, MNRAS, 335, 965 



